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ABSTRACT 

We cross-correlate the third-year WMAP data with galaxy samples extracted from 
the SDSS DR4 (SDSS4) covering 13% of the sky, increasing by a factor of 3.7 the 
volume sampled in previous analysis. The new measurements confirm a positive cross- 
correlation with higher significance (total signal-to- noise of about 4.7). The correlation 
as a function of angular scale is well fitted by the integrated Sachs- Wolfe (ISW) effect 
for LCDM flat FRW models with a cosmological constant. The combined analysis of 
different samples gives Q A = 0.80 - 0.85 (68% Confidence Level, CL) or 0.77 - 0.86 
(95% CL). We find similar best fit values for tl\ for different galaxy samples with 
median redshifts of z ~ 0.3 and z ~ 0.5, indicating that the data scale with redshift as 
predicted by the LCDM cosmology (with equation of state parameter w = —1). This 
agreement is not trivial, but can not yet be used to break the degeneracy constraints 
in the w versus 17a plane using only the ISW data. 



1 INTRODUCTION 



Dark Energy (DE) models with late time cosmic accelera- 
tion, such as the A-dominated CDM model, predict a slow 
down for the growth of the linear gravitational potential at 
moderate redshift z < 1, which can be observed as temper- 
ature anisotropies in the CMB: the so-called late integrated 
Sachs- Wolfe (ISW) effect. The ISW effect is expected to pro- 
duce an increase of power (a bump) in the amplitude of the 
CMB fluctuations at the largest scales, ie lower order multi- 
poles, which are dominated by cosmic variance. This expec- 
tation, seems challenged by observations, as the first year 
WMAP results (WMAP1) confirmed the low amplitude of 
the CMB quadrupole first measured by COBE (eg Hinshaw 
et al. 1996a). The discrepancy between the observations and 
the ACDM model is particularly evident in the tempera- 
ture angular correlation function W2(9), which shows an al- 
most complete lack of signal on angular scales 8 £ 60 de- 
grees. According to Spergel et al. (2003), the probability of 
finding such a result in a spatially-flat ACDM cosmology is 
about 1.5 x 10" 3 . This was questioned in Gaztafiaga et al. 
(2003) who found, using simulated ACDM WMAP maps, 
a much lower significance (less than 2-sigma) for 1112(6). 
A low significance was also estimated by different studies 
(eg Efstathiou 2003, Olivera-Costa et al. 2003), although a 
discrepancy larger than 3-sigma still remains on both the 
quadrupole-octopole alignment (Tegmark, Oliera-Costa & 
Hamilton 2003, Olivera-Costa et al. 2003) and the WMAP 
observed high value of the temperature-polarization cross- 
correlation on large scales (Dore, Holder & Loeb 2004). 



Given the observed anomalies on the ISW predictions, it 
is of particular interest to check if the ISW effect can be de- 
tected observationally through an independent test, such as 
the cross-correlation of temperature fluctuations with local 
tracers of the gravitational potential (Crittenden & Turok 
1996). A positive cross-correlation between WMAP1 and 
galaxy samples from the Sloan Digital Sky Survey (SDSS) 
was first found by Fosalba, Gaztafiaga & Castander (2003, 
FGC03 from now on) and Scraton et al. (2003). FGC03 used 
the lyr WMAP data (WMAP1) and the SDSS data release 
1 (SDSS1). WMAP1 has also been correlated with the APM 
galaxies (Fosalba & Gaztafiaga 2004), infrared galaxies (Af- 
shordi et al 2004), radio galaxies (Nolta et al 2004), and 
the hard X-ray background (Boughn & Crittenden 2004a; 
Boughn & Crittenden 2004b). The significance of these 
cross-correlations measurements was low (about 2-3 a, see 
Gaztahaga, Manera & Multamaki 2006 for a summary and 
joint analysis), and many scientists are still skeptical of the 
reality of these detections. Here we want to check if these re- 
sults can be confirmed to higher significance using the SDSS 
data release 4 (SDSS4) which covers 3 times the volume of 
SDSS1. At the same time, we will compare the signal of 
the 1st and 3rd year WMAP data (WMAP3) recently made 
public (Hinshaw et al 2006; Spergel et al 2006). With bet- 
ter signal-to-noise and better understanding of foreground 
contamination in WMAP3, it remains to be seen whether 
the low significance signal of the WMAP1-SDSS1 analysis 
can be confirmed with WMAP3-SDSS4, or if on the contrary 
this signal vanishes as systematic and statistical errors are 
reduced. 
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2 THE DATA 

In order to trace the changing gravitational potentials we use 
galaxies selected from the Sloan Digital Sky Survey Data 
Release 4 (Adelman-McCarthy et al 2006), SDSS4 here- 
after, which covers 6670 deg 2 (i.e, 16% of the sky). We have 
selected subsamples with different redshift distributions to 
check the reliability of the detection and to probe the evolu- 
tion of the ISW effect. All subsamples studied contain large 
number of galaxies, between 10 6 -10 7 , depending on the sub- 
sample. We concentrate our analysis on the North Galactic 
Cap SDSS4 Area (~ 5500 deg 2 ), because it contains the 
most contiguous area. We have selected 3 magnitude sub- 
samples with r = 18 - 19, r = 19 - 20 and r = 20 - 21 
and a high redshift Luminous Red Galaxy (LRG; e.g. Eisen- 
stein et al. 2001) color selected subsample (17 < r < 21, 
(r - i) > {g - r)/4 + 0.36, (g - r) > 0.72 * (r - i) + 1.7). Be- 
cause of the smaller volume, the r — 18 — 19 and r — 19 — 20 
subsamples provide low signal-to-noise (S/N<2) in the cross- 
correlation with WMAP, and we therefore center our anal- 
ysis on the two deeper subsamples. The mask used for 
these data avoids pixels with observed holes, trails, bleeding, 
bright stars or seeing greater than 1.8. 

To model the redshift distribution of our samples we 
take a generic form of the type: 

N(z) ~ 4> G {z) ~ (*- z c f exp f-£^L) 3/ \ (1) 

for z > z c and zero otherwise. The N(z) distribution of 
the r = 20 — 21 subsample is quite broad with z c — and 
zq ~ 0.2 which results in a median redshift, z — IAzq — 0.3 
(e.g., Dodelson et al. 2001, Brown et al. 2003). On the other 
hand, the LRG subsample has a narrower redshift distri- 
bution. The first colour cut is perpendicular to the galaxy 
evolutionary tracks in the (g-r) .vs.(r-i) colour space and 
ensures that very few z < 0.40 galaxies are selected, which 
translates into a cut z c ~ 0.37 in the above N(z) model. The 
second colour cut is parallel to evolution and perpendicular 
to spectral type differences and selects only red galaxies with 
old stellar populations. The faint magnitude limit (r < 21) 



cuts high redshift galaxies (zo ~ 0.45), which results in an 
overall median redshift of z ~ 0.5. 

We use the full-sky CMB maps from the third-year 
WMAP data (Hinshaw et al 2006; Spergel et al 2006) 
(WMAP3 from now on). In particular, we have chosen the 
V-band (~ 61 GHz) for our analysis since it has a lower 
pixel noise than the highest frequency W-band (~ 94 GHz), 
while it has sufficient high spatial resolution (21') to map 
the typical Abell cluster radius at the mean SDSS depth. 
We use a combined SDSS+WMAP mask that includes the 
KpO mask, which cuts 21.4% of WMAP sky pixels (Bennett 
et al 2003b), to make sure Galactic emission does not af- 
fect our analysis. WMAP and SDSS data are digitized into 
7' pixels using the HEALPix tessellation 1 . Figure shows 
how the WMAP3 and SDSS4 pixel maps look like when den- 
sity and temperature fluctuations are smoothed on 0.3 deg 
scale. 



3 CROSS-CORRELATION AND ERRORS 

We define the cross-correlation function as the expectation 
value of density fluctuations So = No/ < No > — 1 and 
temperature anisotropics At = T — To (in fiK) at two po- 
sitions hi and fi2 in the sky: wtg{9) = (At{&i)Sg(p-2)), 
where 9 = \&z — assuming that the distribution is sta- 
tistically isotropic. To estimate wtg(9) from the pixel maps 
we use: 

Ei,, A T(fil) fc(fij) Wi Wj 

wtg(0) = ^ , (2) 

where the sum extends to all pairs i,j separated by 9 ± A#. 
The weights Wi can be used to minimize the variance when 
the pixel noise is not uniform, however this introduces larger 
cosmic variance. Here we follow the WMAP team and use 

1 Some of the results in this paper have been 
derived using HEALPix (Go rski et al 1998), 
http: / /www. eso.org/science/healpix 
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Figure 2. The continuous line with errorbars shows the 
WMAP3-SDSS4 angular cross-correlation as a function of scale 
for the r = 20 — 21 sample (top) and the LRG sample (bot- 
tom). The dotted line corresponds to using the 1st yr WMAP 
(WMAP1-SDSS4) data, which is very close to the WMAP3 re- 
sults (continuous line). The dashed lines show the ACDM model 
with Q\ = 0.83 ( best overall fit) scaled to the appropriate bias 
and projected to each sample redshift. 
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Table 1. w TG (B) for WMAP3-SDSS4. 



correlation signal. All three estimates give very similar re- 
sults for covariance and the x errors, details will be pre- 
sented elsewhere (Fosalba et al. 2006). 
To compare models we use a \ 2 test: 



JV 

E 

i,j=l 



(3) 



where Ai = w§ G (8i) — u>tg(#i) is the difference between 
the "estimation" E and the model M. We perform a Sin- 
gular Value Decomposition (SVD) of the covariance matrix 



C i 



([/ife) t DkiVij where Dij 



XiSij is a diagonal matrix 
with the singular values on the diagonal, and U and V axe or- 
thogonal matrices that span the range and nullspace of dj . 
We can choose the number of eigenvectors wtg(*) (or prin- 
cipal components) we wish to include in our \ 2 by effectively 
setting the corresponding inverses of the small singular val- 
ues to zero. In practice, we work only with the subspace of 
"dominant modes" which have a significant "signal-to-noise" 
(S/N). The S/N of each eigenmode, labeled by i, is: 



uniform weights (i.e. Wi = 1). The resulting correlation is 
displayed in Fig|2] On scales up to 10 degrees we find signif- 
icant correlation above the estimated error-bars. The dotted 
and continuous lines correspond to WMAP1 and WMAP3 
data respectively, and show little difference within the er- 
rors. This indicates that the cross-correlation is signal dom- 
inated. 

We have used different prescriptions to estimate the co- 
variance matrix: a) jack-knife, b) 2000 montecarlo simula- 
tions c) theoretical estimation (including cross-correlation 
signal) both in configuration and harmonic space. Our mon- 
tecarlo simulations in b) include independent simulations of 
both the CMB and galaxy maps, with the adequate cross- 



( S\ _ WTG(i) _|1 V-rr W TG {j) ... 

\NJi~ Xi -\\ l 2^ U]% a w (j) ■ W 

As S/N depends strongly on the assumed cosmological 
model, we use the direct measurements of wtg to estimate 
this quantity. The total S/N can be obtained by adding the 
individual modes in quadrature. In our analysis we have used 
5 eigenmodes for the r — 20 — 21 sample and 3 for the LRG 
sample. The results are similar if we use less eigenmodes. 
With more eigenmodes, the inversion becomes unstable be- 
cause we include eigenvalues which are very close to zero 
and are dominated by noise. 
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3.1 Comparison with Predictions 

ISW temperature anisotropies are given by (Sachs & Wolfe 
1967): 



A isw f ^ _ T(n) ~T f d$ 

T ^ ' = To = ~ / lz^ ' 



(5) 



where $ is the Newtonian gravitational potential at red- 
shift z. One way to detect the ISW effect is to cross- 
correlate temperature fluctuations with galaxy density fluc- 
tuations projected in the sky (Crittenden & Turok 1996). 
It is useful to expand the cross-correlation w^™ (6) = 
< Ay SW (hi)5a(n2) > in a Legendre polynomial basis. On 
large linear scales and small angular separations it is: 



ISW i a \ 

w TG (6) 



' 1 dzW IS w(z)W G (z)^P(k) 



(21 + 1) 2 
3Ckn(Ho/c 



2 d[D(z)/a] 



dz 



(6) 



Wisw(z) 

W G {z) = b(z)4> G (z)D(z), 

where k — l+1 J 2 , 4>g(z) is the survey galaxy selection func- 
tion in Eq.0 and r(z) is the comoving distance. This is 
just a Legendre decomposition of the equations presented 
in Fosalba & Gaztanaga (2004), see also Afshordi (2004). 
The advantage of this formulation is that we can here set 
the monopole (I = 0) and dipole (I = 1) contribution 
to zero, as it is done in the WMAP maps. The contribu- 
tion of the monopole and dipole to wtg is significant and 
over predicts wtg by about 10%. The power spectrum is 
P(k) = A fc' ls T 2 (k), where T(k) is the ACDM transfer 
function, which we evaluate using the fitting formula of Ein- 
seintein & Hu 1998. 

We make the assumption that on very large scales 
the galaxy distribution is a tracer of the underlaying mat- 
ter fluctuations, related through the linear bias factor, 
Sc(n,z) — b(z)5m(n > z). We estimate b(z) from the angu- 
lar galaxy-galaxy auto-correlation wgg (9) in each sample 
by fitting to the linear flat ACDM model prediction wgg(Q) 
and marginalizing over the value of Q m . The models have 
h = 0.71, Tcmb = 2.725, Q B = 0.022//i 2 , n 3 = 0.938 and 
Qfc = 0. and are normalized to the value of as that best fits 
WMAP3 data (Spergel et al. 2006): <r 8 = 0.75±8;SI- With 
this procedure we find a normalization of bag, ~ 0.90 — 0.96 
and bas ~ 1.02 - 1.12 for the r = 20 - 21 and LRG samples 
respectively. We marginalize all our results over the uncer- 
tainties in both as and bag- This also roughly accounts for 
the uncertainty in the selection function. The predictions 
of wtg do no change much with the selection function (see 
§4.1 in Gaztanaga, Manera & Multamaki 2006), but the bias 
estimated from from wgg depends strongly on the effective 
volume covered by c/>g(z). Because of the marginalization 
our final results do not change much when we change the 
median redshift of the sample by ~ 10%, which represents 
current uncertainties in N(z). But in the case of the LRG it 
is critical to include not only the correct value of the mean 
redshift (or zo in Eq.0) but also the redshift cut z c in- 
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Figure 3. Probability distribution: 1 — P x [> Ax 2 , V = 1] for Ca 
in the r = 20 — 21 sample (short-dashed line), the LRG sample 
(long-dashed line) and the combined analysis (continuous mid- 
dle curve). The range of 68% and 95% confidence regions in 
are defined by the intersection with the corresponding horizontal 
lines. 



troduced by the color selection in Eq.0. In previous LRG 
cross-correlation analysis (eg FGC03 & Gaztanaga, Man- 
era & Multamaki 2006) the value of z c was neglected. This 
can over predict bas, as estimated from wgg, by a factor of 
two. Uncertainties in the shape of N(z) considered here are 
within the normalization errors we have already included for 
as and bas- We have also made predictions for the best fit 
WMAP3 data with n s = 1 which gives different parameters 
and normalization (as = 0.79±o og) and find very similar 
results. 

Under the above assumptions we are left with only one 
free parameter, which is Q m or = 1— fl m . Figl^lshows the 
probability distribution estimated for S7a from the A\ 2 = 
X 2 — Xmin analysis away from the minimum value Xmin ■ Both 
samples prefer the same value of Qa- This is a consistency 
check for the ACDM model. The combined best fit model 
has Q A ^ 0.831^3- The predictions for this S7a best value 
are shown as a dashed line in Fig. [5] 

FigEJshows the joint 2D contours for dark energy mod- 
els with an effective equation of state w = p/p, assuming 
not perturbations in DE and a Hubble equation: H 2 /Hq = 



n(i+zy+n A (i + z 



,3(1+™) 



For each (w, S7a) we derive bas 



consistently from the galaxy-galaxy auto-correlation data. 
We also marginalize over the uncertainties in bas and over 
as 6 (0.65,0.85), to account for the WMAP3 as normaliza- 
tion for to/-l. The cosmological constant model w = — 1, 
however, still remains a very good fit to the data. This is due 
to the large degeneracy of the equation of state parameter w 
with Oa- This degeneracy can be broken by supernovae SNIa 
data (eg see Corasaniti, Giannantonio and Melchiorri 2005 
and Fig. 8 in see Gaztanaga, Manera & Multamaki 2006). 
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Figure 4. Two dimensional contours for Ha and w, the DE 
effective equation of state. The inner black contour limits the 
ID marginalized 68% confidence region (Ax 2 = 1). The other 
contour correspond to 95% limits (Ax 2 = 4). 



4 DISCUSSION 

The objective of our analysis was primarily to check if we 
could confirm or refute with higher significance the findings 
of the WMAP1-SDSS1 cross-correlation by FGC03. With 
an increase in area of a factor of ~ 3.7 in SDSS4, larger 
signal-to-noise and better understanding of foregrounds in 
WMAP3, our new analysis shows that the signal is robust. 
This is also in line with the first findings using optical 
(APM) galaxies by Fosalba & Gaztanaga (2004). The cross- 
correlation signal in WMAP3-SDSS4 seems slightly larger 
than in previous WMAP1-SDSS1 measurements which re- 
sults in slightly larger values for £Ia (see FGC03 and 
Gaztanaga, Manera & Multamaki 2006). This is probably 
due to sampling variance, as the SDSS4 volume has increase 
by almost a factor of 4 over SDSS1. We find little difference 
within the errors in the cross-correlation of WMAP1-SDSS4 
and WMAP3-SDSS4 (see FigHJ. 

The total S/N in Eqg|of the WMAP3-SDSS4 correla- 
tion is S/N ~ 3.6 for the r = 20 - 21 sample and S/N ~ 3.0 
for the LRG, which gives a combined S/N ~ 4.7, assum- 
ing the two samples are independent. We have checked the 
validity of this assumption by doing a proper join analysis 
where we include the covariance between the two samples. 
For the join analysis we find a S/N ~ 4.4 with the first 2 
dominant eigenvectors and S/N ~ 4.8 with 4 eigenvectors. 

We find that a ACDM model with Qa — 0.83 success- 
fully explains the ISW effect for both samples of galaxies 
without need of any further modeling. The best fit for Oa 
for each individual sample are very close. This is signifi- 
cant and can be understood as a consistency test for the 
ACDM model. 

The equation of state parameter appears to be very 
degenerate and it is not well constrained by current ISW 
data alone (see also Corasaniti, Giannantonio and Mel- 
chiorri 2005 and Gaztanaga, Manera & Multamaki 2006). 
Upcoming surveys such as the Dark Energy Survey (DES, 
www.darkenergysurvey.org), with deeper galaxy samples 



and more accurate redshift information should be able to 
break the w — J1a degeneracy and maybe shed new light on 
the the nature of dark energy. 
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